KA-TP-18-2001 



tiep-ph/0105349 



The Implementation of the Minimal Supersymmetric Standard Model 

in FeynArts and FormCalc 

Thomas Hahn, Christian Schappacher 
Institut fur Theoretische Physik 

Universitdt Karlsruhe 
D-76128 Karlsruhe, Germany 

June 1, 2001 
Abstract 

We describe the implementation of the MSSM in the diagram generator FeynArts and 
the calculational tool FormCalc. This extension allows to perform loop calculations of MSSM 
processes almost fully automatically. The actual implementation has two aspects: The MSSM 
Feynman rules are specified in a new model file for FeynArts. The computation of the param- 
eters in the MSSM Lagrangian from the input parameters is realized as a Fortran subroutine 
in the framework of FormCalc. The model file does not depend on the latter, however, and can 
be used even if one does not want to continue the calculation with FormCalc. The Feynman 
rules have been entered in a very generic way to allow e.g. scenarios with complex param- 
eters, and have been tested extensively by reproducing known results for several non-trivial 
scattering processes. 

PACS numbers: 12.60.Jv, 02.70.-c. 

1 Introduction 

One of the main problems of Feynman-diagrammatic computations is the enormous growth of 
the number of Feynman diagrams, not only with the loop order, but also with the number of 
particles in a model. While many precision calculations in the Standard Model (SM) could still be 
performed by hand exactly, the same is very difficult in models like the Minimal Supersymmetric 
Standard Model (MSSM). Yet it is highly desirable to perform unabridged calculations in the 
MSSM, too, since also the MSSM allows to make precise predictions in terms of a set of input 
parameters. 

With the availability of powerful software packages, the basic problem of bookkeeping and 
calculation of the diagrams has been solved for many common cases. Still, it is not entirely trivial 
to code a model of the complexity of the MSSM in such a system, since this has to be done in a 
reasonably general way (i.e. not only for special cases of the parameters) and many checks have 
to be performed to test all sectors of the model. 

The present paper documents the implementation of the MSSM in the FeynArts and Form- 
Calc fj| packages. Other programs for which an MSSM model file exists are GRACE || |^| and 
CompHEP (|, 0]. Conceptually, the FeynArts-FormCalc system works in three stages, as sketched 
in the following diagram. 
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A more detailed discussion of the interplay between FeynArts and FormCalc can be found in |7j]. 

This paper is organized as follows. Corresponding to the three stages, Diagram generation, 
Algebraic simplification, and Numerical evaluation, there are three sections which describe the 
modifications to each stage that are necessary for calculations in the MSSM. 

Sect. ^] describes the changes in FeynArts, namely the new MSSM model file. The model 
file declares the properties of the fields, their propagators, and their couplings. It contains the 
parameters of the MSSM Lagrangian. The Feynman rules have been entered in a very generic way 
to allow e.g. scenarios with complex parameters. 

Sects. H and [| describe the changes in FormCalc. Sect. |^ outlines the changes in the algebraic 
simplification and Sect. |] explains the calculation of the parameters that are used in the model 
file. The model-file parameters are derived from a (reasonably small) set of input parameters. 
This is realized as a model-initializing Fortran subroutine in the framework of FormCalc. The 
model file does not depend on the latter, however, and can be used even if one does not want to 
continue the calculation with FormCalc. 



2 Diagram generation 

2.1 Algorithms 

Supersymmetric theories contain fermion-number-violating couplings, e.g. quark-squark-gluino, 
so that the usual method of ordering the Dirac matrices oppositely to their occurrence along the 
arrows on fermionic lines breaks down since one cannot define a fermion-number flow. 

FeynArts uses the "flipping-rule" algorithm || . This algorithm was invented precisely to solve 
the problem of fermion-number-violating couplings and works as follows: Instead of traversing 
the fermion lines along the fermion-number flow imposed from the outside, FeynArts chooses a 
direction for each fermion chain. If it turns out later that, for a Dirac fermion, the chosen direction 
is opposite to the actual fermion flow, FeynArts "flips" the coupling, i.e. it derives the coupling 
appropriate for the reversed fermion flow from the known coupling. The flipping of a coupling is 
in fact nothing but a charge (as opposed to hermitian) conjugation. 

2.2 The MSSM model file 

The model file is the source of all physics information in FeynArts. It declares the properties of the 
fields, their propagators, and their couplings. In the model file the parameters of the Lagrangian 
are used, not a restricted set of input parameters. 

There are two versions of the MSSM model file in FeynArts, both of which follow the conven- 
tions of [§, [L0[ |ll|. The file MSSMQCD . mod defines the complete (electroweak and strong) MSSM, 
whereas MSSM. mod contains only the electroweak subset, defined as everything except the gluon, 
its ghost, and the gluino. The four-sfermion couplings appear in MSSM. mod although they have 
both electroweak and strong parts. Counter-terms are not yet entered in the model files. 

Table [j] gives the names for the fields and corresponding masses defined in MSSM. mod and 
MSSMQCD. mod. The symbols used for the MSSM parameters are specified in Table |^. The complete 
list of couplings is too long to be included here, but is contained in the FeynArts distribution as 
a PostScript file. 

The MSSM model files further declare a number of restrictions with which certain groups of 
particles can be excluded in the diagram generation. These restrictions are listed in Table |. 



3 Algebraic simplification 

3.1 Algorithms 

The algebraic simplification has to address two problems that arise in one-loop calculations in 
supersymmetric theories in general and in the MSSM in particular. The first is the conceptual 



2 



leptons 


f = P 


field 


mass 


sleptons 


f = P 


field 


mass 


v a 




F[l,{g>] 





"a 




S[ll, {ff}] 


MSf 


tg 




F[2,{g}] 


MLE 


c g 




S[12, {s, 5 }] 


MSf 



quarks 


squarks 


Ug 




F[3, {g,o}] 


MQU 






S[13, {s, 5 ,o}] 


MSf 


dg 




F[4,{<?,o}] 


MQD 






S[14, {s,3,o}] 


MSf 



gauge bosons 


neutralinos, charginos 


7 


yes 


V[l] 





Xn 


yes 


F[ll,{n}] 


MNeu 


Z 


yes 


V[2] 


MZ 


Xc 




F[12, {c}] 


MCha 


w- 




V[3] 


MW 











Higgs bosons 


ghosts 


h° 


yes 


S[l] 


MhO 


u~ ( 




U[l] 





H° 


yes 


S[2] 


MHH 


uz 




U[2] 


MZ 


A 


yes 


S[3] 


MAO 


u + 




U[3] 


MW 


G° 


yes 


S[4] 


MGO 


U- 




U[4] 


MW 


H~ 




S[5] 


MHp 






U[5, {u}] 





G~ 




S[6] 


MGp 











gluon 


gluino 


.9 


yes 


V [5 , {«}] 





9 


yes 


F[15, {u}] 


MG1 



where the following indices are used: 



g = Index [Generation] = 1 . . . 3 , s = Index [Sf ermion] = 1 . . . 2 , 

o = Index [Colour] = 1 . . . 3 , n = Index [Neutralino] = 1 . . . 4 , 

u = Index [Gluon] = 1 . . . 8 , c = Index [Chargino] = 1 . . . 2 . 



Table 1: The particle set-up in MSSM.mod and MSSMQCD .mod. The gluon, its ghost, and the gluino, 
which are defined only in the latter, are written in grey. 
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MhO, MHH, MAO , MGO 
MHp, MGp 
CB, SB, TB 
CA, SA 

C2A, S2A, C2B , S2B 
CAB, SAB, CBA, SBA 
MUE 


neutral Higgs masses 
charged Higgs masses 
cos/3, sin/3, tan/3 
cos a, sin a 

cos 2a, sin 2a, cos 2/3, sin 2/3 

cos(a + p), sm(a + p), cos(/3 — a), sm(p — a) 

Higgs-doublet mixing parameter \i 


MG1 

MNeu [n] 
ZNeu[n, n'] 
MCha [c] 

UClia[c, c'] , VCha[c, c'] 


gluino mass 
neutralino masses 
neutralino mixing matrix 
chargino masses 
chargino mixing matrices 


MSf [s, t, 

USf [i.fif] [s, s'] 

Af [t, 5 ] 


sfermion masses 

sfermion mixing matrix 

(scalar) soft-breaking A-parameters 



The indices enumerate the following properties: 
sfermion generation: g = 1 ... 3 , 



sfermion type: t 



1 sneutrinos , 

2 sleptons , 

3 up-type squarks , 

4 down-type squarks , 



eigenstates: mass 

sfermion s = 1 . . . 2 , 

neutralino n = 1 . . . 4 , 

chargino c = 1 . . . 2 , 



gauge 
a' = 1 ... 2 , 
ri = 1 . . . 4 , 
c' = 1 . . . 2 . 



Table 2: Symbols representing the MSSM parameters introduced by MSSM.mod and MSSMQCD.mod. 



NoGenerationl 
NoGeneration2 
NoGeneration3 


exclude generation- 1 fcrmions (V , e, u, d) 
exclude generation-2 fcrmions (z/^, /*, c, s) 
exclude generation-3 fermions (v T , r, t, b) 


NoElectronHCoupling 
NoLightFHCoupling 


exclude all couplings involving electrons and any Higgs field 
exclude all couplings between light fermions (all fermions 
except the top) and any Higgs field 


NoSUSYParticles 
THDMParticles 


exclude the particles not present in the SM: sfcrmions, 
charginos, neutralinos, and the Higgs fields H°, A , 

exclude the particles not present in the two-Higgs-doublet 
model: the sfermions, charginos, and neutralinos 



Table 3: Pre-defined restrictions in MSSM.mod and MSSMQCD.mod. 
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problem of implementing a supersymmetry-preserving regularization scheme, the second is the 
technical problem of dealing with more diagrams and larger expressions. 

The default regularization scheme employed by FormCalc is dimensional regularization. Un- 
fortunately, this scheme is known to break supersymmetry FormCalc has therefore been 
equipped with an alternative scheme suited for calculations in supersymmetric theories. It is the 
constrained differential renormalization scheme (CDR) which is equivalent to regularization by 
dimensional reduction at the one-loop level ||. The schemes are chosen with the Dimension option 
of the CalcFeynAmp function: Dimension -> D selects dimensional regularization (the default), 
Dimension -> 4 switches to CDR. 

The second problem is a technical one. In the SM, both the number of diagrams and the 
comparatively simple coupling structure allow to perform calculations with "brute force." To wit, 
the diagrams are generated with all indices explicitly written out, and then calculated. 

The situation is less favourable in the MSSM, however: not only is the number of diagrams 
considerably higher in most cases, but also the coupling structure is more involved. This is because 
particles like sfermions, or gauginos and higgsinos in general mix to form mass eigenstates and 
hence their couplings contain plenty of mixing matrices and can become rather lengthy. This 
means that the algebraic simplification has to be performed more carefully in order to maintain a 
decent performance and keep the size of the results as small as possible. 

The way in which FormCalc proceeds is already suggested by the level structure of the Feyn- 
Arts amplitudes. The lowest (generic) level completely determines the kinematical structure of 
a diagram, so FormCalc first performs the time-consuming simplifications involving kinematical 
quantities, like the tensor reduction, on the far fewer generic diagrams only. For each diagram at 
the next higher (classes) level, it then substitutes the generic coupling constants and masses by 
their actual values. The diagrams are mostly specified by then, except that index summations, 
e.g. over fermion generations, are not yet carried out. These index summations are performed only 
in the Fortran code, which means that FormCalc has to keep track of all indices and write out 
the corresponding loop instructions in Fortran. For a model like the MSSM, the savings incurred 
with this method of simplification may easily amount to an order of magnitude in both CPU time 
and size of the Fortran code. 

4 Numerical evaluation 

In the FormCalc framework, the Mathematica expressions resulting from the algebraic simplifi- 
cation are translated to Fortran code for numerical evaluation. The generated Fortran code of 
course has to be provided with the proper numerical values for the parameters appearing in the 
model, i.e. the variables in Table |^. This is solved by a subroutine which is called at the beginning 
of the calculation to initialize all model parameters. 

Technically, the input parameters are specified in common blocks defined in model . h. Some less 
commonly changed inputs, notably some breaking parameters in the sfermion sector, are realized as 
preprocessor variables which take default values if not defined by the user. Preprocessor variables 
are also used for switches, such as whether squark mixing should be turned on or off. 

4.1 Parameters of the MSSM 

Supersymmetry (SUSY) completely determines the supersymmetric part of the MSSM Lagrangian 
once the SM parameters are known. The complete MSSM, however, with softly broken SUSY in 
its general form, necessarily introduces a large number of masses, phases, and mixing angles to 
parametrize the SUSY breaking. At a final count, 105 of these degrees of freedom cannot be 
absorbed in some other quantity or rotated away |L4|. In order to get a handle on so many 
parameters, several assumptions have been made in mssm_ini.F to arrive at a moderate number 
of input parameters while retaining reasonable generality. 

The first assumption is that the SUSY-GUT relation holds, which relates the U(l), SU(2), and 
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SU(3) gaugino mass parameters Mi, M 2 , and M 3 according to 

Mi = l^M 2 and M 3 = ^V s w M 2 with m§ = |M 3 | , (1) 
3(% a(s) 

where yfs is the CMS energy, cw = Mw/Mz, and sw = \/ 1 — c w ■ ^ ne running couplings a and 
a s are MS values; we use the CERNlib function ALPHAS2 JD| for a s (s) and 

1 1 J_V^ n 2 w1 _J_ 1 20 s 

a(s) ~ a(M 2 z ) ~ Sir ^ Qf c Mf ~ 5(M§) ~ 9tt M| ( j 

with l/a(M|) = 127.934 |§. 

Secondly, we assume the SUSY-breaking parameters in the sfermion sector to be flavour-blind 
and distinguish only two v4-parameters for different isospin values, hence the sfermion sector is 
governed by seven scalar parameters: 



(3) 



where in addition the various M^- get the default value MgusYj a common SUSY-breaking mass, 
i.e. 

Mq =M l = My =M b = M g = A/susy (by default) . (4) 
The remaining MSSM input parameters thus have 10 degrees of freedom: 

1 quotient of VEv's: tan f3 (real), 

1 Higgs mass: Ma (real), 

6 breaking parameters: A u ,Ad (complex), MsuSY,-^ (real), 

2 mixing parameter: \x (complex). 

The inputs to mssm_ini . F are summarized in Table [I| All other variables in the MSSM Lagrangian 
(see Table ||) are determined once this reduced set of MSSM input parameters and the SM inputs 
are specified and are calculated by mssm_ini . F. This is the reason why mssm_ini . F is much more 
involved than its companion file for the SM, sm_ini .F. The calculation of the parameters can be 
directed by several switches, which are also listed in Table ||. 

When the particle masses are calculated from these input parameters, mssm_ini.F checks 
whether they are consistent with the current exclusion limits and automatically omits already- 
excluded points with a warning. The following bounds are used: 



(5) 



m~ t Ss 80 GeV (17) , m > f 91 GeV for real input Parameters |1 



Ss 70 GeV g7| , ' h ' |85 GeV for complex input parameters |§ , 

™&b,i > 150 GcV . 90 GeV , (6) 



m l > 70 GeV || , m . > 30 GeV || , 
A Pl - b < 3 x 10~ 3 @ , m~ g > 175 GeV || . 

The detailed formulas built into mssm_ini.F for computing the various parameters are pre- 
sented in the following sections, organized according to the various sectors of the MSSM. 

4.2 The Higgs sector 

The neutral Higgs sector is fixed by choosing a value for tan/3 = v-ijvi (the ratio of the vacuum 
expectation values of the two Higgs doublets) and for the mass Ma of the CP-odd neutral Higgs 
boson A . For the CP -even Higgs masses, which receive sizable radiative corrections, we use the 
approximation formula of |23) which agrees with the full two-loop calculation 0] to within less 
than 2 GeV. 
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parameter 


Fortran name 


tVDC 


default value 


tan f3 


TB 


double precision 




M A 


MAO 


double precision 




A u 


Au 


double complex 




A d 


Ad 


double complex 






MSusy 


double precision 




M 2 


M_2 


double precision 






MUE 


double complex 






MSQ 


preprocessor variable 


MSusy 


Mi 


MSL 


preprocessor variable 


MSusy 


M 


MSU 


preprocessor variable 


MSusy 


M 3 


MSD 


preprocessor variable 


MSusy 


M E 


MSE 


preprocessor variable 


MSusy 



switch 


action 


NO_SQUARK_MIXING 


Sets A u — fi*cotp and Ad = ^*tan(3, so that the off-diagonal 
entries of the sfermion mass matrices vanish (cf. Eq. (|37|)), i.e. 
makes sfermion mass eigenstates = sfermion gauge eigenstates. 


COMPLEX_PARAMETERS 


Uses a simpler (one-loop) approximation for the Higgs masses 
which is valid for all parameters. This switch must be set if com- 
plex input parameters are used because otherwise the more precise 
(two-loop) approximation for the Higgs masses is taken, which is 
valid only for real parameters. 


SM_DNLY 


Calculates Mh as usual, but then reverse-engineers the mixing in 
the Higgs sector (a and (3) such that the MSSM Higgs sector looks 
like a SM Higgs sector (see 11 , p. 356), with the light CP-even 
Higgs boson h figuring as the SM Higgs boson. 


NO_EXCLUSION_LIMITS 


Ignores the experimental bounds, i.e. does not exclude points in 
parameter space if the bounds in Eq. (|^) are violated. 



Table 4: The inputs and switches for mssm_ini . F. Note: A preprocessor variable is defined with 
"#define var value" and a switch is set with "#define switch". The #define must stand at the 
beginning of the line. 
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4.2.1 The neutral CP-even Higgs bosons h and H 

In this section we discuss the case of real MSSM parameters, for which a more sophisticated 
computation of the Higgs masses is implemented. The case of complex parameters is treated in 
Sect. |4~2~2] 

The gauge eigenstates 4>\ and fa of the neutral CP-even Higgs bosons mix via the mass matrix 

M 2 = ( M\ sin 2 /? + Ml cos 2 /3 - E 01 -(M 2 + Mf) sin/3 cos /3 - E 010 2 \ 
Higgs \-(M% + M%) sin $ cos l3-% l4 , 2 M\ cos 2 /3 + Mf sin 2 /3 - E 02 y 

Diagonalizing this matrix yields the Higgs masses 



2 

1 



^i + M^ + (E 01 ^E^_ MlM|cos22/3 



4 



-,1/2 

+ | (Ml - M 2 )cos2/3(E 01 - E 02 ) + (M^ + Mf) sin2/3E 0102 + E 2 



1>1<32 



(8) 



The mixing angle a of the CP-even Higgs doublet is hence 



a = arctan ( + M 2 ) sin/^ - E^ 2 \ 

\ M 2 Z cos 2 /? + M 2 sin 2 /? - E^ -Ml) 

The presence of renormalized self-energies E indicates that the mass matrix contains significant 
radiative corrections which must be taken into account to make quantitatively correct predictions. 
The self-energies of course receive contributions from all sectors of the MSSM, but not all are 
numerically of equal relevance. We take into account only the following terms: 

• Yp-J/t) (0): the one-loop i/i-contributions at zero momentum transfer up to mf p5{ , 

• E( 2 '*/^(0): the dominant two-loop i/i-contributions of 0(aa s ) at zero momentum transfer 
|23f and the leading two- loop Yukawa correction of 0(a 2 ) |f26f , 

• E( 1,rest ): the one- loop leading-log contributions from all other sectors p7| . 

In the following, we give the explicit expressions for these contributions as implemented in 
mssm_ini . F. 



The one-loop i/i-contributions at zero momentum transfer up to m 
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(13) 
(14) 
(15) 



• The dominant two- loop i/i-contributions of 0(aa s ) at zero momentum transfer and the 
leading two-loop Yukawa correction of 0(a 2 ): 
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(18) 



and {/ is the stop mixing matrix defined in Eq. 

rn t = m t (m t ) 



We use the MS top mass 



l + £a s (m 2 ) 



(19) 



instead of the pole mass to include also the leading t/t-contributions beyond 0(aa s ). 
• The one-loop leading-log contributions from all other sectors: 
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v„ 



4.2.2 The case of complex parameters 

The sophisticated calculation of the Higgs masses as in the preceding section is valid only for real 
parameters. When complex parameters are used, the preprocessor variable COMPLEX_PARAMETERS 
must be defined (see Table ^), which makes mssm_ini .F use a simple (one-loop) approximation 
for the Higgs masses that is valid also for complex parameters. 

For the Higgs masses the dominant one-loop contributions are the m 4 terms p5|. Neglecting 



the splitting of the t masses, the 02 self-energy of Eq. ( p"2| ) becomes 



s(l,t/t) 



(0) 
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(24) 



with the running top mass rat as defined in Eq. (|19|). All other radiative corrections in the Higgs 
mass matrix are set to zero. 

4.2.3 The charged Higgs bosons 

There are only small radiative corrections for the charged Higgs masses. We use the following 
relation, valid for M A ~ 0(M Z ) |1, up to M^-values of 250 GeV: 
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(25) 



with N r . 



N„ 



3 and the running top mass m t as defined in Eq. (|19|). Above 250 GeV, this 



formula reduces to M~ ± — M\ 



M 2 ,. 



4.3 The chargino sector 

The charged gauginos and higgsinos mix via the mass matrix ||, ^] 

M 2 V2M w s\n(3\ 
\/2Mw cos P [i I 



X = 



(26) 
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X is diagonalized with two unitary matrices U and V according to 

[/*Xyt = diag ( m - i)m - 2 ) ) (27) 

which yields the chargino masses 



Y 



(29) 



< 2 = Ml±l^±?Mk T j M + W + MbY _ {M , w sin w _ ^ 2 . (28) 
For the numerical diagonalization, the LAPACK subroutine ZGESVD [§8] is used. 

4.4 The neutralino sector 

The neutral gauginos and higgsinos mix via the mass matrix |l0) 

/ Mi -Mzs w cos/3 M z s w smP\ 

M 2 M z cw cos /3 -M z c w sin /3 

—Mzsw cos /3 Mzcw cos /3 — /i 

\ M z siy sin /3 —Mzcw sin /3 -/i / 

V" is diagonalized with a unitary matrix A according to 

A* y = diag(m^.o , m^o , m^o , m^o ) , (30) 

which yields the four neutralino mass eigenstates. For the numerical diagonalization, the LAPACK 
subroutine ZGESVD p3| is used. Note that we work with a general complex matrix N and also 
check that the computed mass values are nonnegative. 

4.5 The sfermion sector 

The sfermion gauge eigenstates are connected via the mass matrix || |l0[ 



Mf L + mj m f (Mf R )* 
m f Mf- R Mf R + ml 



(31) 



where 



Z 3 W 1^1 f° r left-handed sleptons, 

M ifl = A - ^* {° 0t/3 f ° r U_typ6 sfermions ( J 3 = + 1 / 2 ), / 33 s 
' [tan/3 for <i-type sfermions (7g = —1/2), 

{M^ for right-handed, u-type squarks, 
M? for right-handed, d-type squarks, (34) 
M| for right-handed sleptons. 

The mass eigenstates are obtained by diagonalizing Z with a unitary matrix {/* , viz. 

U f ZU" = diag(m| , mjj , (35) 

which results in the sfermion masses 

to 2 = to 2 - 1 - — ' l ll-T^R i . / ( AjLL 



„„2 + ± j^ M LL + M fl« -p ^ ( M ll _ m «r^ + Am)\Mf R \ z j . (36) 

For the numerical diagonalization the LAPACK subroutine ZHEEV |2^] is used. 

Sfermion mixing can be switched off by defining the preprocessor variable NO_SQUARK_MIXING 
(see Table ^) . This causes the breaking parameters A u and Ad to be set to the values 

A u = /J* cot j3 , A d = fi* tan (3 (37) 

so that the off-diagonal terms in the mass matrix Z vanish. 
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4.5.1 Supersymmetric contributions to Ap 



In addition to the experimental bounds on the squark masses, mssm_ini . F also checks whether 
the MSSM contributions to the p-parameter are consistent with the current exclusion limits (see 
Eq. @). This bound becomes relevant mainly when parameters are chosen to achieve a large 
mass splitting between t\ and ti or b\ and hi. 

We have implemented the calculation of the MSSM contributions to Ap according to [ p9| . 
The sfcrmion contributions are significant only when the masses of the isospin partners are very 
different, and since this is governed by the quark masses (see Eq. (|3^)), only corrections from 
t/6-loops have been included. Contributions from gluino exchange, which are very lengthy and 
vanish for large m g , have similarly been neglected. 

These contributions to Ap are given by 

Apf,g = UtiUl*UilU&F{ml,ml) + ^^^C/^ (m^m^) 

+ WtM^Fim^ml) + mM^Fimlml), (38) 
where the are the sfermion mixing matrices defined in Eq. ( |35| ) and the function 

*-, r , 3Gp _ . . Gp ce s (ui}) _ . . ,„„% 

F ^ y) = ^ F ^ y) + w^* F2ix > y) (39) 

includes the one- and two-loop contributions 

F l (x,y)=x + y-^-ln-, (40) 
x-y y 

F 2 (x, y =x + y y - 2 + - In - In - + §■ In 2 - - 2 (x - y Li 2 1 - - . 41 

x-y \ y yj y (x-yy y \ y J 

The two-loop contribution is of the order of 10 to 15% of the one-loop result. 



5 Tests 

The model files MSSM. mod and MSSMQCD.mod have been checked against results known from the 
literature for a variety of scattering processes. In cases where the Fortran programs of the original 
authors were available, we found perfect agreement for all differential cross-sections. In the other 
cases we could reproduce qualitatively the figures. The list of processes we checked together with 
the agreement we achieved is given in Table §. 



6 Availability, Requirements 



The FeynArts package including the MSSM. mod and MSSMQCD.mod model files can be downloaded 
from http://www.feynarts.de. The program itself requires Mathematica 3 or above. For the 
topology editor, which is not necessarily invoked, a Java VM and the J/Link package are needed, 
both of which can be obtained free of charge (see the instructions on the web site) . 

The FormCalc package is available from |http : //www. f eynarts .de/f ormcalc and includes 
the initialization file mssm_ini.F. It runs on Unix-like platforms and requires Mathematica 3, 
FORM 3, and the GNU C compiler (gcc). To link the Fortran code generated by FormCalc, one 



needs in addition the LoopTools library (tittp : //www. f eynarts . de/looptools) and the CERNlib 



(tittp : //wwwinf o . cern. ch/ asd/index .html 



FeynArts and FormCalc each include a comprehensive manual which explains installation and 
usage. Both are open-source programs and stand under the GNU library general public license. 
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process 


agreement 


reference 


scope 


qq — > tt 


11 digits 




SUSY-QCD only 


qq — > tt 


10 digits 




full MSSM 


55 ^ 


qualitative 


m 


SUSY-QCD only 


55 W 


10 digits 


J In 

m 


full MSSM 


5ff -> H+H- 


10 digits 


n 


full MSSM 




qualitative 


J In 

H 


full MSSM 


e+e- -> 


11 digits 


n 


full MSSM 


e+e" -> W+W~ 


qualitative 


Hi 


sfermion contributions 


e+e- — ^ ff+fl"- 


7 digits 


@ 


full MSSM 


dipole moments 


11 digits 




full MSSM 



Table 5: Checks performed with FeynArts and FormCalc to test the implementation of the MSSM. 
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